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Abstract. —Sequence variation was examined in the cytochrome b region of the mitochondrial genome 
of Burrowing Owls {Athene cunicularia) from North and South America, and compared with the Elf Owl 
{Micrathene whitneyi). Barred Owl {Strix varia) and Eastern Screech Owl {Otus asio). Attempts to clone 
and sequence the control region of the mitochondrial genome resulted in sequences that appeared to 
be nuclear copies of that region. Cytochrome b sequences revealed a genetic split between Burrowing 
Owl populations from North and South America. This split may date back 2 million yr to the connection 
of these continents via the isthmian land bridge. Additional population structure appears to be of 
Pleistocene origin or more recent. Data indicate a possible North American origin for Burrowing Owls 
and subsequent dispersal via the land bridge to the South American continent. The depth of the split 
between Burrowing Owls from North and South America is consistent with species-level distinction. 
Additional data from nuclear markers, morphology and/or ecological indicators, such as behavior or 
vocalizations, will be necessary to confirm these results. 

Key Words: Burrowing Owl; Athene cunicularia; mitochondrial DNA; cytochrome t, genetics; North America; 
South America. 


Analisis preliminar de la estructura del adn mitocondrial en poblaciones de Buho Cavador 

Resumrn. —La variacion de la secuencia fue examinada en la region del citocromo b del genoma mi¬ 
tocondrial de Buhos Cavadores {Athene cunicularia) de Norte y Sur America, y fue cornparada con las 
del Buho Elfo {Micrathene whitneyi), cl Buho Barreteado {Strix varia) y el Buho Chirriador oriental {Otus 
asio). Los intentos para clonar y secuenciar la region de control del genoma mitocondrial dieron como 
resultado secuencias que parecian ser copias nucleares de esa region. Las secuencias del citocromo b 
revclaron una division genetica entre las poblaciones de Buhos Cavadores de Norte y Sur America. Esta 
escision puede datar de 2 millones de ahos atras cuando se conectaron estos dos continentes por medio 
del puente terrestre del istmo. Otra estructura de la poblacion parece tencr origen en cl Pleistoceno o 
mas recientemente. Los datos indican un posible origen Norteamericano para cl Buho Cavador y una 
subsecucnie dispersion hacia el continentc Suramericano a traves del istmo. La ptofimdidad de la 
division entre los Buhos Cavadores de Norte y Sur America es consistente con el nivel de distincion a 
especie. Datos adicionales a partir de marcadores nucleares, morfologia y/o indicadores ecologicos, al 
igual que comportamiento o vocalizacioncs, scran necesarios para conlirmar estos resulLidos. 

[Traduccion de Victor Vanegas y Cesar Marquez] 

southern Canada through Argentina, including 
the West Indies. Suitable habitat includes arid ar¬ 
eas, savannas, and grasslands (Pregill and Olson 
1981, Johnsgard 1988). Burrowing Owls nest in 
burrows and, depending on geographic locality, 
nest in burrow systems of colonial sciurids, other 


The Burrowing Owl {Athene cunicularia) is widely 
distributed throughout the New World, from 
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burrowing vertebrates, natural cavities, or dig their 
own burrows. Currently treated as a single species, 
the Burrowing Owl has 18 recognized subspecies 
described on the basis of plumage characteristics 
and geographic variation in size (Peters 1940). 

The Burrowing Owl fossil record, although 
sparse, and the biogeographic history of savanna/ 
arid habitats throughout the New World have led 
to the development of hypotheses about the evo¬ 
lution of this species. A presumed owl ancestor 
Speotyto megalopeza first appeared in the fossil re¬ 
cord in the Pliocene in western Meade County, 
Kansas, U.S.A. (Ford 1966). The first recorded Bur¬ 
rowing Owl fossil from South America is from the 
Pleistocene (Vuilleumier 1985). The most exten¬ 
sive fossil record of the Burrowing Owl is found in 
the West Indies from the Pleistocene (Pregill and 
Olson 1981). The savanna/arid habitats evolved 
separately on the North and South American con¬ 
tinents during the Cenozoic. For most of this pe¬ 
riod, the two continents were physically separated 
and the centers for the evolution for arid land 
communities were located centrally within each 
continent (Webb 1977). In North America, the 
amount of open habitat replacing forests increased 
substantially throughout the Cenozoic. Today, 
greater than 25% of the land cover of North Amer¬ 
ica is nonforest (Webb 1977). The full extent of 
savanna and arid habitats in South America during 
this same period is not fully understood and re¬ 
mains controversial. Evidence indicates the devel¬ 
opment of flora associated with open habitats dur¬ 
ing the early Cenozoic. The presence of open 
country vertebrates from fossil sites dates to the 
same period and includes early ungulates, chin- 
childs, octodontids, and large grassland birds, such 
as rheas {Rhea spp.) and the carnivorous phorus- 
rhacids (Webb 1978). Climatic shifts of alternating 
humid and dry periods in South America, associ¬ 
ated with glacial and interglacial periods, contrib¬ 
uted to the disjunct distribution of savannas, grass¬ 
lands, and other xeric habitats (Haffer 1974). This 
likely contributed to the disjunct populations of 
Burrowing Owls currently found throughout the 
South American continent. 

Based on the current disjunct distribution of 
some owl populations, climatic fluctuations, and 
evidence of historical extinctions and range reduc¬ 
tions of Burrowing Owls in the West Indies, we hy¬ 
pothesized that mitochondrial DNA (mtDNA) 
would reveal population structure by broad geo¬ 
graphic locality. Specifically, we predicted in North 


America the resident Florida subspecies {A. c. flor- 
idana) would be genetically distinct from western 
populations {A. c. hypugaea) because of its geo¬ 
graphic isolation. We also predicted a genetic 
break would exist between populations found in 
North America and South America. Within South 
America, we anticipated finding genetic differenc¬ 
es between broad areas because of the apparent 
restricted range and isolation of many populations 
and their presumed resident status. 

Methods 

To test our hypotheses, samples from five recognized 
subspecies of Burrowing Owls were analyzed. Blood sam¬ 
ples of <0.5 cc were collected from the western subspe¬ 
cies {A. c. hypugaea) in Nebraska, South Dakota, Oregon, 
California, and New Mexico. Blood samples of the Flor¬ 
ida subspecies (A. c. floridana) were collected by Brian 
Mealy (Univ. Florida-Miami). Tissue samples were ob¬ 
tained from museum tissue collections for Burrowing 
Owls from Baja California, Mexico (A. c. hypugaea, N — 
1); Providences Tucaman and Corrientes, Argentina (A 
c. cunicularia, N = 2)‘, Providence Trujuilo, Peru (A c 
nanodeSy N = 1); and Providence Loja, Ecuador (A 6 
punensis, N = 1). The initial analysis for broad geographic 
variation was conducted on one specimen from each of 
the widely separated populations from Nebraska, Califor¬ 
nia, New Mexico, Florida, Baja California, Mexico, Ar¬ 
gentina (Provinces Tucaman and Corrientes), Peru, and 
Ecuador to evaluate the potential of the selected marker 
to identify population structure. Samples were also ob¬ 
tained from the Elf Owl {Micrathene whitneyi) , Barred Owl 
{Strix varia), and Eastern Screech Owl {Otus asio) for 
comparative purposes. 

Target DNA Sequences. Our original objective was to 
examine two regions of the mitochondrial genome (the 
control region and a section of the cytochrome b region) 
to address geographic variation within subspecies, as well 
as broad-scale variation. However, attempts to use south¬ 
ern blotting to develop specific control region primers 
to address finer-scale questions resulted in sequences that 
appear to be nuclear copies of this region (Desmond 
1997). These sequences are therefore not discussed fur¬ 
ther in this paper. Universal cytochrome b primers de¬ 
signed by Kocher et al. (1989) were used to amplify the 
297 bp section from the cytochrome b region. Results 
from these sequences are the focus of this paper. 

DNA Extraction. DNA was extracted from blood and 
tis.sue by digestion in 100 mM Tris pH 7.5, 100 mM EDTA 
pH 8.0, 100 mM NaCl, 0.5% SDS, and Proteinase K (0.5 
ug/ml) overnight at 55°C. The DNA was purified by ex¬ 
tracting once with phenol, chloroform (24:1, chloro- 
form:isamyl alcohol) and twice with chloroformiisamyl al¬ 
cohol (24:1). Each sample was then ethanol precipitated 
overnight, pelleted and washed twice with 70% ethanol, 
dried and re-suspended in a buffer of IM Tris and lOM 
EDTA. 

Polymerase Chain Reaction (PCR) Amplification. Am¬ 
plification was performed in 50 ul reactions containing 
10 niM Tris-HCl (pH 8.3), 50 niM KCl, 1.5 mM MgClg, 
each dNTP at ImM, each primer at 5 uM, and 2 units of 
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Table 1. Numbers and percentages (in parentheses) of the four nucleotides in the cytochrome b sequences of the 
Burrowing Owl, Barred Owl, and Elf Owl, respectively. 



A 

C 

T 

G 

Athene cunicularia 

71 (24) 

106 (36) 

70 (24) 

49 (17) 

Micrathene whitneyi 

82 (28) 

96 (33) 

71 (24) 

43 (15) 

Stnx varia 

76 (26) 

99 (33) 

71 (24) 

-50 (17) 

Otus asio 

73 (25) 

102 (35) 

70 (24) 

51 (17) 


Taq. Each cycle of the polymerase chain reaction con¬ 
sisted of denaturation at 94°C for 1 min, annealing at 
47°C for 1 min, and extension for 1 min 30 sec at 72°C. 
This cycle was repeated 35 times. All cytochrome b am¬ 
plifications were sequenced using an automated sequenc¬ 
er at the Iowa State University DNA sequencing facility. 
Sequence results were checked manually using the pro¬ 
gram Editview (Version 1.0, Perkin Elmer). 

To test against the possibility of contamination, blood 
samples of Burrowing Owls from Nebraska, Florida, Cal¬ 
ifornia, and New Mexico were taken to a separate labo¬ 
ratory, where no avian work was being conducted. Using 
all new stock solutions, DNA was purified and amplified. 
PCR reactions were conducted under UV light as an ad¬ 
ditional precaution against contamination. All unique se¬ 
quences were deposited in Genbank under the following 
accession numbers: 

nkit 423326, 423334, 424445, 424887, 423676, 424455, 

423688, 423690, 424905, 424907, 423692. 

Results 

Cytochrome b Sequences. A single 297 bp frag¬ 
ment of the cytochrome b gene was consistently 
amplified for Burrowing, Elf, Barred, and Eastern 
Screech Owls, all of which were sequenced. Using 
the same primers, fragments of similar size have 
been amplified for other bird species (Kocher et 
al. 1989, Edwards and Wilson 1990, Birt-Friesen et 
al. 1992). Compared to sequences deposited in the 
Genbank, owl sequences were most similar to other 
avian cytochrome b sequences. The nucleotide 
content of the cytochrome b sequences were con¬ 
sistent with other published avian cytochrome b se¬ 
quences showing a lower guanine (G) content and 


above average cytocine (C) content (Table 1). Also 
in agreement with other avian cytochrome b se¬ 
quences, a third codon deficiency of G and thy¬ 
mine (T) was observed (Table 2). 

Maximum parsimony analysis of the cytochrome 
b data was conducted with PAUP (version 3.1.1; 
Swofford and Begle 1993) using the heuristic 
search algorithm (Elf, Barred, and Eastern Screech 
Owls were designated as outgroups). Unweighted 
maximum parsimony analysis separated Burrowing 
Owls from North and South America, and boot¬ 
strap analysis (500 replications) provided strong 
support for this separation with a value of 100%. 

Cytochrome B Intraspecific Sequence Diver¬ 
gence. Within the 297 bp segment of the cyto¬ 
chrome b region, 1.5 positions were variable among 
all samples (Appendix). With one exception, vari¬ 
able positions were transitions. There were 11 sub¬ 
stitutions in third codon positions between North 
and South American Burrowing Owls. There was 
one substitution in the first position of a codon 
between North American Burrowing Owls (A. c. hy- 
pugaea and A. c. floridana) and A. c. cunicularia 
from Argentina, and two substitutions were in the 
first position of a codon between North American 
Burrowing Owls and A. c. nanodes and A. c. punen- 
sis, which resulted in one amino acid change (Ap¬ 
pendix) . The western (A. c. hypugaea) and Florida 
(A. c. floridana) Burrowing Owls each differed 
from the Peru (A. c. nanodes) and Ecuador {A. c 


Table 2. Percent nucleotide distribution at first, second, and third codon positions of 297 bp cytochrome b fragment 
for Burrowing, Elf, Barred, and Eastern Screech Owls. 




First 



Second 



Third 


A 

G 

C 

T 

A 

G 

C 

T 

A 

G 

C 

T 

Athene cunicularia 

23 

26 

24 

27 

16 

20 

27 

34 

26 

7 

59 

7 

Micrathene whitneyi 

26 

24 

26 

25 

20 

19 

25 

37 

38 

3 

48 

10 

Stnx varia 

24 

25 

24 

26 

24 

18 

28 

36 

28 

7 

55 

9 

Otus asio 

23 

26 

28 

23 

21 

19 

24 

36 

28 

7 

53 

9 
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Table 3. Pairwise distances (based on cytochrome b data) among taxa. Above the diagonal are absolute distances 
and below the diagonal are total number of observed differences. Total number of observed differences equals 
number of transitions plus transversions; number of transversions are in parentheses. 


Taxa 

1 

2 

3 

4 

5 

6 

7 

8 

A, c, hypugaea 

— 

0.007 

0.040 

0.044 

0.044 

0138 

0.185 

0.178 

2. A. a floridana 

2 (0) 

— 

0.040 

0.044 

0.044 

0.138 

0.178 

0.178 

3. A, c. cunicularia 

12 (0) 

12(0) 

— 

0.010 

0.010 

0.1.35 

0.168 

0.168 

4. A. c. nariodes 

13(1) 

13(1) 

3(0) 

— 

0.000 

0.135 

0.168 

0.168 

5. A. c. punensis 

13(1) 

13 (1) 

3(0) 

0(0) 

— 

0.135 

0.168 

0.168 

6. Micraihene whitneyi 

41 (19) 

41 (19) 

40 (19) 

40 (19) 

40 (19) 

— 

0.135 

0.152 

7. Otus asio 

55 (22) 

55 (22) 

50 (23) 

50 (23) 

50 (23) 

40 (16) 

— 

0.135 

8. Strix varia 

53 (19) 

53 (19) 

50 (16) 

50 (16) 

51 (16) 

45 (19) 

40 (18) 

— 


punensis) forms by 12 transitions and one transver¬ 
sion (percent divergence: P = 4.4) and the Argen¬ 
tinean form {A. c. cunicularia) by 12 transitions {P 
= 4.0). The western and Florida Burrowing Owls 
differed by two transitions (percent nucleotide di¬ 
vergence: P = 0.7). The Argentinean Burrowing 

Western North America 



Figure 1. Geographic localities of the distinct mtDNA 
genotypes: A = Athene cunicularia hypugaea, B = A. c. flor- 
idana, C = A. c. nanodes and A. c. punensis, D = A. c. 
cunicularia. 


Owls (A. c. cunicularia) differed from Burrowing 
Owls of coastal Ecuador (A. c. punensis) and Peru 
(A. c. nanodes) by three transitions {P = 1.0), and 
no differences were observed between A. c. nanodes 
and A. c. punensis (Table 3). 

Rates of Evolution. Wood and Krajewski (1996) 
estimated the rate of evolution of cytochrome b 
sequences for cranes (Sarus spp.) to be 1.7% per 
1 million yr, lending support to the molecular 
clock calibration of 2% sequence divergence per 
million yr by Brown et al. (1982). Applying this 
conventional mammalian mtDNA molecular clock 
calibration of 2% sequence divergence per million 
yr to the cytochrome b sequence data, we estimate 
that the North and South American forms of the 
Burrowing Owl diverged ca. 2 million yr ago. Of 
the three South American subspecies examined, A 
c. cunicularia diverged from A. c. punensis and A. c 
nanodes ca. 500 000 yr ago, and coastal populations 
representing A. c. nanodes and A. c. punensis of Peru 
and Ecuador showed no evidence of divergence. 
The two North American subspecies (A. c. hypugaea 
and A. c. floridana) diverged more recently, ca. 
350000 yr ago (Eig. 1). This molecular clock has 
been widely applied to various taxa; however, given 
uncertainties about mtDNA clocks in birds and 
other taxa (Rising and Avise 1993), this estimate 
should be cautiously interpreted. Rates of evolu¬ 
tion were not determined among owl species using 
cytochrome b data due to the saturation of nucle¬ 
otide substitutions. 

Discussion 

Cytochrome b Sequences. The cytochrome b se¬ 
quences obtained for Burrowing Owls and other 
owl species were in agreement with other pub¬ 
lished avian cytochrome b sequences (Kocher et al. 
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1989, Birt-Friesen et al. 1992, Wood and Krajewski 
1996). Like other avian cytochrome b sequences, 
all four owl species exhibited a deficiency of gua¬ 
nine. When nucleotide distribution was examined 
in relation to codon position, the four owl se¬ 
quences exhibited a pattern similar to other avian 
taxa, exhibiting a deficiency for both guanine and 
thymine at third codon positions. The guanine de¬ 
ficiency also has been reported in rodents and fish 
(Kocher et al. 1989). A deficiency in both guanine 
and thymine at the third codon position has only 
been reported for avian species (Kocher et al. 
1989, Quinn and Wilson 1993). 

The level of variability detected within and be¬ 
tween Burrowing Owl populations from North and 
South America, using cytochrome b, is higher 
than, but consistent with, the level of cytochrome 
b variability detected in other avian studies that 
make intraspecific comparisons (Edwards and Wil¬ 
son 1990, Wenink et al. 1993). The agreement in 
observed nucleotide composition with other avian 
species, and the similarity in levels of observed var¬ 
iability within other avian species suggests that the 
cytochrome b sequence data for the four owl spe¬ 
cies is authentic mitochondrial cytochrome b se¬ 
quence. In addition, the amino acid translations 
for the nucleotides did not reveal any stop codons 
interrupting translation that would indicate a non¬ 
functional copy. 

Distribution and Divergence of Cytochrome b 
Sequences. Cytochrome b sequence results re¬ 
vealed a clear split between Burrowing Owls from 
North and South America (4.0-4.4% divergence). 
The sequence data contained numerous diagnostic 
sites for North and South American Burrowing 
Owls. The depth of this split is consistent with spe¬ 
cies-level distinction that has been described in 
other studies using this molecular marker (Ed¬ 
wards and Wilson 1990, Birt-Friesen et al. 1992, 
Wood and Krajewski 1996) and studies using re- 
sliiclion endonucleases on intDNA (Kessler and 
Avise 1984, Mack et al. 1986, Shields and Wilson 

1987) . Within-genus sequence divergence among 
Saras cranes (Wood and Kajewski 1996) is similar 
to divergence values between North and South 
American populations of Burrowing Owls. Several 
studies examining sibling taxa (Kessler and Avise 
1984, Shields and Wilson 1987, Avise and Zink 

1988) report lower values for genetic distances 
than observed between North and South American 
forms of the Burrowing Owl. The one notable ex¬ 
ception is the study on Australian Babblers {Po- 


matostomus temporalis Edwards and Wilson 1990), 
which reports large within-subspecies values for ge¬ 
netic distance using cytochrome b. They report a 
mean of 1.4% and 1.3% within northern and 
southern groups of P. temporalis, respectively, and a 
mean of 3.2% between the northern and southern 
forms. They comment that this degree of diver¬ 
gence is large for subspecific status and suggest fur¬ 
ther investigation into the specific status of the 
bird. The even greater mtDNA distances between 
the North and South American forms of the Bur¬ 
rowing Owl reported here suggest that the species 
status of these owls should also be reevaluated. The 
distances observed between A. c. hypugaea and A. 
c. floridana for North America and between A. c. 
cunicularia and A. c. nanodes/punensis for South 
America are consistent with subspecies-level dis¬ 
tinctions that have been observed in other studies. 

Although a strong split between Burrowing Owls 
from North and South America was detected, we 
cannot geographically define exactly where the 
split occurs because of our limited sampling. It 
most likely occurs in Central America or extreme 
northern South America (i.e., Colombia). A more 
intensive sampling effort between the areas in 
question would determine whether or not this 
break is between two distinct clades. 

Biogeographic Patterns. The divergence of 
North and South American forms of the Burrow¬ 
ing Owl ca. 2 million yr ago coincides with the 
presence of the isthmian land bridge, which arose 
at that time, providing a dispersal corridor be¬ 
tween continents. The mingling of the North and 
South American faunas is postulated to have 
spanned a 2 million yr period between the last mil¬ 
lion years of the Pliocene and the first million years 
of the Pleistocene (from 3 until 1 million yr ago). 
Webb (1978) estimated that a large portion (2/3 
for mammals) of the faunal exchange involved sa¬ 
vanna-adapted organisms. The presence of the 
Burrowing Owl ancestor in the North American 
fossil record dating back to the Pliocene, and 
South American fossil records of the Burrowing 
Owl from the Pleistocene, may indicate that these 
birds evolved in North America and subsequently 
dispersed to the South American continent; how¬ 
ever, the sparse fossil record may be misleading. 
The climatic fluctuations in South America associ¬ 
ated with glacial and interglacial periods resulted 
in isolated patches of suitable habitat that ranged 
from small to large expanses of grasslands and sa¬ 
vanna-type habitats. The disjunct distribution of 
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suitable habitat throughout South America helps 
explain the diversity among Burrowing Owl popu¬ 
lations on that continent (11 recognized subspe¬ 
cies). The genetic divergence observed between cy¬ 
tochrome b sequences of Burrowing Owls from 
North and South America is large and warrants fur¬ 
ther investigation regarding full species status. Ver¬ 
ification of the observed cytochrome b sequence 
divergence with nuclear markers, morphology 
and/or other ecological indicators, such as vocali¬ 
zations or behavior, is needed. 
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Appendix. Nucleotide and inferred amino acid sequences of 297 bp fragment of the cytochrome b gene for Bur¬ 
rowing Owl subspecies and the Elf, Eastern Screech, and Barred Owls. Sequence orientation is from 5' to 3' on the 
light strand. Dots indicate identity with A. c. hypugaea sequences, and dashes indicate gaps in sequences. Sites in bold 
are amino acid replacements within Burrowing Owls, and sites in italics indicate amino acid replacements among 
owl species. 
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Appendix, Continued. 
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